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1. Introduction 

The pseudospin-electron model (PEM) was originally proposed to include a local 
interaction of the conducting electrons in metals (or semimetals) with a certain two 
level subsystem when it is reasonable to use pseudospin formalism; the pseudospin 
variable Sf = ±1/2 defines these two states. Starting from |l|, |^ this scheme is apphed 
to describe the strongly correlated electrons of Cu02 sheets coupled with the vibrational 
states of apex oxygen ions Oiv (moving in a double-well potential) in YBaCuO type 
systems. Recently a similar model was proposed for investigation of the proton-electron 
interaction in molecular and crystalline systems with hydrogen bonds [P, Q . The model 
Hamiltonian has the following form: 

i ijcr 

and includes, in addition to the terms describing electron transfer (~ tij), the electron 
correlation [U - term) in the spirit of the Hubbard model, interaction with the 
anharmonic mode {g - term), the energy of the tunnelling splitting {Q - term), and 
the energy of the anharmonic potential asymmetry {h - term) in the single- site part 

Hi = Un^riii - /i(ni| + rijj + g{n^ + rijJS'i' - VtS^ - hS-. (2) 

Based on PEM a possible connection between the superconductivity and lattice 
instability of the ferroelectric type in HTSC was discussed P, ^. Description of the 
electron spectrum and electron statistics of the PEM was given in within the 
framework of the temperature two-time Green's function method in the Hubbard-I 
approximation. 

A series of works has been carried out in which the pseudospin (S^S^), mixed (S^n), 
and charge (nn) correlation functions were calculated. A possibility of divergences of 
these functions at some values of temperature exists as it was shown with making use of 
the generalized random phase approximation (GRPA) |^ in the limit of infinite single- 
site electron correlations {U —>■ oo) p, ^. This effect was interpreted as a manifestation 
of dielectric instability or ferroelectric type anomaly. The tendency to form the spatially 
modulated charge and pseudospin ordering at certain model parameter values was 
estabhshed. 

The case of absence of electron transfer {tij = 0), but with taking into account 
of the direct interaction between pseudospins was considered within the mean field 
approximation [[Toi [ll], It was shown there that behaviour of the system strongly 
depends on a thermodynamical regime: at yU = const the first order phase transitions 
with the jumps of pseudospin mean value {S^) and electron concentration values n as 
well as the second ones (with smooth changes) were detected, while at n = const an 
instability with respect to phase separation in the electron and pseudospin subsystems 
can take place. 



In the papers [0, |Tj, |15|, |1^ PEM was considered at the absence of direct 



electron-electron interaction and tunnelling splitting {U = 0, Q = 0). In such a 
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version of PEM (simplified PEM), an effective many-body pseudospin interaction via 
conducting electrons appears, and hence it is interesting to investigate an influence of the 
retarded nondirect interaction between pseudospins on thermodynamics of the model. 
Hamiltonian of this simplified PEM is invariant with respect to the transformation 
/i — >■ — /i, /i 2(7 — /i, n — i> 2 — n, S*^ — > —S^ (the so-called electron-hole symmetry) 
what makes possible to describe the hole-pseudospin system as well. 

On the other hand, the presented Hamiltonian in the case of absence of electron 
correlation and tunnelling splitting allows one to describe the binary alloy type model. It 
is convenient to introduce projective operators on pseudospin states = l/2±S'f ; then 
Hamiltonian of the binary alloy can be obtained by substitution P^ = pi, P[^ = 1 — 
where Pi is the concentration of one component of binary alloy, and 1 — pi is the 
concentration of the second one. Difference between these models is in the way how 
an averaging procedure over projection operators P^ is performed (thermal statistical 
averaging in the case of PEM and configurational averaging for binary alloy) and how 
the self-consistency is achieved (fixed value of longitudinal field h for PEM and fixed 
value of the component concentration p for binary alloy). 

Furthermore, if in Hamiltonian (|I]), (^ (considering the simplified case U = 0, 
f2 = 0), we remove spin indices and rewrite the Hamiltonian in terms of the operators 
of the mobile di {ci„ = di) and localized fi {P^ = f^fi, P~ = 1 — f^fi) electrons, 
we obtain the Hamiltonian of the Falicov-Kimball model where h plays a role of the 
chemical potential for the localized /-electrons. The ground state of the Falicov-Kimball 
model, when the electron concentration for subsystems is fixed, is not uniform and shows 
the commensurate or incommensurate ordering or phase separation depending on the 



concentration and coupling constant values [|T^]. On the other hand, in the case of the 
fixed value of field h, the bistability effects are observed. 

Investigations of the thermodynamic properties of the simplified PEM within the 
framework of the self-consistent GRPA [T^ and dynamical mean field approximation 
(DMFA) have shown that: 

• in the /i = const regime (when the electron states of other structure elements, 
which are not included explicitly into the PEM, are supposed to play a role of 
a thermostat ensuring a fixed value of the chemical potential fi) the interaction 
between the electron and pseudospin subsystems leads to the possibility of either 
first or second order phase transitions between different uniform phases (bistability 



effect) [[131 , |l^ as well as between the uniform and the chess-board ones [[151 , |T6|; 

in the regime n = const (this situation is more customary at the consideration 
of electron systems and means that the chemical potential depends now on the 
electron concentration being the function of T, h etc.) an instability with respect 
to the phase separation in the electron and pseudospin subsystems can take place 



The above mentioned methods of investigation of the PEM (within the framework 
of the self-consistent GRPA scheme as well as within the framework of DMFA in the 
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case of the finite space dimension) take into account only mean field type contributions. 
It is reasonable only when fluctuation effects are small. In the vicinity of critical points 
effects of the mean field fluctuations become significant. 

The aim of this article is to calculate the thermodynamic and correlation functions 
of the simplified PEM ([/ = 0, ^2 = 0) within the framework of self-consistent approach, 
allowing us take into account fluctuation effects of the effective self-consistent mean field. 
Such a consistent gaussian fluctuation approach was proposed previously in work 



where expressions for the grand canonical potential and pseudospin correlator as well as 
a self-consistent set of equations for the pseudospin mean value and root-mean-squares 
(r.m.s.) fluctuation parameter were obtained. 

Furthermore we would like to investigate an influence of the gaussian fluctuations on 
the thermodynamics of phase transitions in the PEM and discuss an applicability of the 
schemes previously used for the investigation of the PEM. For this purpose the results 
of investigation of PEM within the framework of the self-consistent GRPA scheme 



as well as within the DMFA [|14[ are presented. 

The paper is organized as follows. First, a short review of the self-consistent 
GRPA approach is presented and the results of the numerical calculations within the 
framework of this scheme are compared with corresponding ones obtained within the 
DMFA in W^. Second, we review the self-consistent gaussian fluctuation approach for 



the PEM and present its simplified version (analogous to the approximation proposed 
by Onyszkiewicz |20[ for spin models). Then we apply gaussian fluctuation method 
to investigation of the thermodynamic properties of the PEM, and finally results of the 
numerical calculations of the pseudospin mean value and r.m.s. fluctuation parameter 
as well as the grand canonical potential are shown. Phase diagrams are presented and 
compared to the corresponding ones of the GRPA. Finally, discussion and conclusion 
are given. 

2. Self-consistent GRPA method 

We would like to remind that traditional GRPA method (proposed by Izyumov et al. 
for the investigation of the magnetic susceptibility of the ordinary Hubbard model as 
well as t — J model) takes into account the same topological class of diagrams (so-called 
loop diagrams) as in the random phase approximation (RPA). Main difference between 
the GRPA and RPA is a way of choosing of the basic states: splitted (Hubbard-I) band 
states in the GRPA method and pure band states in RPA. The question how to calculate 



thermodynamic quantities within the GRPA has been open until recent works [0, [21 
In this section we briefly show how to construct the mean field type approximation 
within the GRPA scheme for PEM. 

Calculations are performed in the strong coupling case {g ^ t) using single-site 
states as the basic ones. The formalism of electron annihilation (creation) operators 
O'ia = Ci^Pt ^ ciia = CiaPf (P^ = 1/2 ± S^) actlug at a site with the certain pseudospin 
orientation is introduced JTSf. Within the framework of the proposed representation. 
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the initial Hamiltonian (|I]) (in the case of f/ = 0, f2 = 0) has the following form: 

h 
2 



H = Ho+Hir,t = ^{einia+e2nia-^S-} + ^ tijiatraja+atraj^+dl^aj^+dl^dj^), (3) 



where £i_2 = — /i ± (7/2 are the energies of the single-site states. 

Expansion of the calculated quantities in terms of the electron transfer i^int leads 
to the infinite series of terms containing the averages of the T-products of the ai„, di^, 
af^, df^ operators. The evaluation of such averages is made using the appropriate Wick's 



theorem [|13l formulated in the spirit of Wick's theorem for Hubbard operators This 
theorem gives an algorithm reducing the average of product of n creation (annihilation) 
operators to the sum of averages of products of the n — 2 operators. So it is possible 
to express the result in terms of the products of nonperturbed Green's functions and 
averages of the projection operators expanded in semi-invariants. 
Nonperturbated electron Green's function is equal to 

P+ pr 

gM = {9iM); 9iM = -. — h : — - — . (4) 

In the uniform case {S-) = {S^), a single-electron Green's function (calculated 
in Hubbard-I type approximation) is = Gk{ujn) = (fi'^^(t^n) — tk)'^', its poles 

determine the electron spectrum 



,(tfc) = i(-2/x + tfc) ± + Uk{S^)g + tl . (5) 



Investigation of the electron spectrum was performed in In this approximation, 

the branches Siitk) and Suitk) form two electron subbands always separated by a gap. 
Reconstruction of the electron spectrum takes place with change of the pseudospin mean 
value {S^)- 

In the adopted approximation the diagrammatic series for the pseudospin mean 



value can be presented in the form [13 




(5^) = ©-(:6:^^+^,c^:^^^^-... , (6) 

where the following diagrammatic notations are used: Q ~ "S*^, — C::=0— ~ Qiif^nj^ wavy 
line is the Fourier transform of the hopping integral t^. Semi-invariants are represented 
by ovals and contain the (5-symbols on site indices. In the spirit of the traditional mean 
field approach, renormalization of the basic semi-invariant by insertion of independent 
loop fragments is taken into account in (^). The analytical expression for the loop is 
the following: 

N - I - 




^ =P{arPt + a2Pn- (7) 



It should be noted that within the presented self-consistent scheme of the 
GRPA, chain fragments form the single-electron Green's function in the Hubbard-I 
approximation, and in the sequences of loop diagrams in expressions for thermodynamic 
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and pair correlation functions the connections between any two loops by more than one 
semi-invariant are omitted. This procedure includes renormalization of the higher order 
semi-invariants similar to the one given by expression (^). 

Summation of diagrammatic series (^) for the pseudospin mean value is equivalent 
to averaging with the mean-field type Hamiltonian: 



Hmf= ^{ei{ni^+nii) + 62(711^ +nii) - yS-}, 



(8) 



where an expression y = h + a2 ~ ai determines an internal effective self-consistent 
field formed by retarded many-body interaction between pseudospins via conducting 
electrons. 

Equation for pseudospin mean value in the uniform case is as follows 

(5f) = (5>F = %), (9) 

where 



%) = -tanh|-|/ + ln^^-^ 



(10) 



An analytical expression for mean value of the particle number can be presented in the 
form [HI: 



(11) 



= ]7 E + n{eM)\ - [1 + 2{S^)]n{e2) - [1 - 2{S^)]n{e,). 



Here nie] 



1 



is Fermi distribution. 



1 + e^^ 

Diagrammatic equation on pseudospin correlator {S^S^)q within the framework of 
the GRPA has the following form: 




(12) 



Equation ([T2|) differs from the corresponding one for the Ising model in RPA by the 
replacement of the direct exchange interaction by the electron loop Ilq = 

(describing the many-body retarded interaction between pseudospins via conducting 
electrons) 



n,k 



in{k) 



;i - Qntk) 



(13) 



The first term in equation (|T2]) is the zero-order correlator renormalized due to the 
inclusion in all basic semi-invariants of the mean-field type contributions ('single-tail' 
parts) coming from the effective pseudospin interaction 



(14) 




and is thus calculated with the help of Hamiltonian i^MF- 
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In analytical form, solution of equation (|1^) is equal to 
1/4- (^^)^ 
1 + j^Y.KUk)tn{k + q)(l/4 - {S^Y 



n,k 

and is different from zero only in a static case {un = 0) (since pseudospin operator 
commutes with the initial Hamiltonian) . 



Within the same approach the grand canonical potential was obtained ||13[. The 
corresponding analytical expression is 

OCKP. ^ -^5:in(l-*..K)) - ^5:^|=M^^ - ^,„T.(e-»»). (16, 

n.k n,k 

It was shown that all quantities can be derived from the grand canonical potential 



d(-/i) ' " d{-h) ' " d{(3h) 
what demonstrates thermodynamic consistence of the proposed approximation. 

In the /i = const regime the stable states are determined from the minimum 
of the grand canonical potential. The solution of equation for the pseudospin mean 
value and calculation of potential Q (as well as the pair correlation function) were 
performed numerically [jl3|. The first order phase transitions between the different 
uniform phases (bistability effect) at change of temperature T and field h can take 
place when chemical potential is placed within the electron subbands. In the case 
when chemical potential is fixed and placed between the electron subbands, the uniform 
phase become unstable with respect to fluctuations with q = (tt, tt), and the possibility 
of second order phase transitions between the uniform and chess-board phase exists at 
the change of temperature or field ||16 |. 

As mentioned above, the band structure is determined by the pseudospin mean 
value (figure P, and its change is accompanied by the corresponding changes of the 
electron concentration (see |jl3|, |15| for details). 

The simplified version of PEM was considered in [jl^ within the framework of 
DMFA scheme as well. The obtained phase diagrams within the self-consistent GRPA 
and DMFA are presented in figure |^ (in the case when chemical potential is placed in 
the center of lower subband /i = —0.5). 

One can see a quite sufficient coincidence of shapes of the first order phase transition 
curves (thick solid line). But unlike to the phase diagram in figure ^ obtained within 
the DMFA, boundaries of the phase stability region calculated in the self-consistent 
GRPA have the same type of slope, hence a region exists, where the vertical line twice 
crosses the boundary of the phase stability region (figure 0a). The analysis of the 
(S^S^) behaviour in this region (for fixed value of the longitudinal field h = 0.0875) 
with the decrease of temperature shows that high temperature phase becomes unstable 
with respect to fluctuations with q = (the lower crossing point of the vertical line 
and boundary of the phase stability region in figure Ha.) Similar results were obtained 




Figure 1. Field dependences of pseudospin mean value and electron bands boundaries 
for fj, = const regime when chemical potential is placed in the lower subband (in the 
self-consistent GRPA) {g = 1, tk=o = 0.2, ^ -0.5, T = 0.0132) 
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Figure 2. — h phase diagrams: (a) within the self-consistent GRPA, (b) within 
the DMFA 10. Solid and dashed lines indicate the first order phase transition line 
and boundaries of the phase stability (spinodal lines), respectively (g — 1, tk=a = 0.2, 
fi = -0.5) 



previously in 0, |^ for temperature behaviour of the correlation functions in the case of 
infinite single-site electron correlation U ^ oo within the framework of the GRPA. 

In the case of the fixed value of the electron concentration condition of equilibrium 
is determined by the minimum of free energy. In this regime the first order phase 
transition with a jump of the pseudospin mean value (accompanied by the change 
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of electron concentration) transforms into a phase separation into the regions with 
different electron concentrations and pseudospin mean values |jT^. For the first time 
the instability with respect to the phase separation in pseudospin-electron model was 
marked in [Q, where it was obtained within the GRPA in the limit U 



oo. 



3. Self-consistent gaussian fluctuation approximation 

The above considered approach takes into account only contributions of mean field 
type. In this section we present the developed in [|l8l consistent scheme for calculation 
of thermodynamic and correlation functions, which takes into account the gaussian 
fluctuations of the self-consistent mean field. We also reduce the presented here method 
to more suitable for numerical calculation taking into account a restricted class of 
diagrams. Such an approximation for spin models with a direct interaction proposed 
by Onyszkiewicz |l^, yields a much better description of critical properties of spin 
models in the whole range of temperature then others. 

In constructing a higher order approximation that takes into account the fluctuation 
effects of the self-consistent mean field, we use the self- consistent GRPA as the zero-order 
one. This means that all 'single-tail' parts of diagrams @ and ([T^) are already summed 
up and all semi-invariants are calculated using the distribution with the Hamiltonian 
Hmf (§)• AH these semi-invariants are represented graphically by thick ovals and contain 
the 5-symbols on site index: 



--{S^)mf = Ky), 



c 
MF 



dPy ' 



/ QZ QZ QZ\C 



[17] 



As an approximation that goes beyond of the self-consistent GRPA we use, similarly 
to [^, the approach taking into account the so-called 'double-tail' diagrams. A 
corresponding diagrammatic series for the pseudospin mean value can be written as 




The diagrammatic equation for the pseudospin correlator {S^S^)q within the presented 
here approximation is similar to the corresponding one within the GRPA (JT^) and given 

by 




(19) 



but now all semi-invariants in this equation are renormalized due to the 'double-tail' 
parts. Thus the corresponding diagrammatic series on zero-order correlator looks like 
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The contribution corresponding to the 'double-tail' fragment of the diagrams can be 
written in the following analytical form (using the notation (p!3|)): 

iV3 



X = f-iri = ^5^5^5^A^t„(fc)t.(fc-q)(S^S^),A2i„,(/c')t„'(fc'+g)- 

n,n' k,k' Q 

The equation on pseudospin correlation function (pj]) has such a solution: 



(21) 



1 + f EA2t„(fc)t„(fc+g)H 

n,k 



(22) 



Since the pseudospin correlator (|2^ ) is a frequency independent function, in the 
expression for X (^1]) we have two independent sums over internal Matsubara 
frequencies, allowing one (by decomposition into simple fractions) to sum over all 
internal frequencies. 

The diagrammatic series (|l^) and (^) can be expressed in the following analytical 



forms [18 



V2^ 



exp 



2X 



(23) 



1 



XV2^ 



exp 



2X 



(24) 



As one can see the contribution of diagrammatic series with 'double-tail' parts 
corresponds to the average with the Gaussian distribution, where X can be interpreted 
as the root-mean-squares (r.m.s.) fluctuation of the mean field around the mean value 
of y. Thus we obtain a self-consistent set of equations for pseudospin mean value (|^) 



and r.m.s. fluctuation parameter (^) as well as expression for pseudospin correlation 
function (|2^). 

The grand canonical potential within the approximation presented here is given by 
the diagrammatic series below: 




11 

21! 

The corresponding analytical expression is 

+ 00 

" — "GRPA + TT 



(25) 



j e-^e%+Ode-^ j {l-erf(-|L)}sign(0&(z/+Ode 



2 V27rX 

— oo — oo 

E ^ninik? + ^ In (i + s| 5^ Aiukr) . 



(26) 



n,k 



n,k 
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The grand canonical potential written in this form satisfies the stationary 
conditions: 

^ = ^ = (27) 



which are equivalent to the equations (^) and (pSl). The consistency of the 
approximations used for the pseudospin mean value, pseudospin correlation function 
and grand canonical potential can be checked explicitly using the relations: 

T? iTT ^ ^'^^)' un\ = {S^S^)q=0. (28) 

d[—h) a{—h(j) x=const 

In the limit of vanishing fiuctuations our results go over into the ones obtained within 
the self-consistent GRPA. 

The analytical scheme presented here for the pseudospin-electron model can be 
easily reduced to the scheme proposed by Onyszkiewicz for spin models. For this purpose 
we consider renormalization with use of the simplest possible pseudospin correlation 
function involving only gaussian fiuctuations of the mean field for the 'double-tail' 
fragment of the diagram 




Within the framework of this simplification the grand canonical potential satisfies the 
stationary conditions (|27|) and can be written in the following analytical form: 

+ 00 

1^ = l^GRPA + ^HX - i J |l-erf(-|L)}sign(0% + Ode (30) 

— oo 

The diagrammatic series for pseudospin mean value is the same as the above presented 



ones (18). Finally, the set of equations on pseudospin mean value and r.m.s. fiuctuation 



parameter can be written down as: 

+ 00 



(S') = -jL= f exp f-l^ ]Hy + Od(, (31) 



4. Phase diagrams within the gaussian fluctuation approach 

Solution of equations for the pseudospin mean value and r.m.s. fiuctuation parameter 
is performed numerically for the square lattice with nearest-neighbor hopping. The 
stable state of the system is described by the solution (from the possible set of ones) 
corresponding to the global minimum of grand canonical potential; metastable states 
are related to solutions corresponding to local minima. 
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Figure 3. T^ — h phase diagrams: (a) within the self-consistent GRPA, (b) within the 
gaussian fluctuation methods, (g — 1, ifc^o — 0.2, ji — —0.4). 



We would like to note that there is no particular difference (figure ^p) between 
the results obtained within the self-consistent gaussian fluctuation approach by solving 
the set of equations (PT|), ( PB| ) and its simplified version (the set of equations (PT|), 
(p^)): the same topological type and slope of phase diagrams, similar field and 
temperature behaviours of pseudospin mean value and grand canonical potential were 
observed. Only small quantitative differences (as one can see in phase diagrams in 
figure are seen (the thin solid phase coexistence line is obtained in the self-consistent 
gaussian fluctuation approach, thick solid line corresponds to the Onyszkiewicz type 
approximation). Therefore, to show the influence of the gaussian fluctuations on the 
thermodynamic properties of the PEM we perform all calculations for the simplified 
variant of the gaussian fluctuation approach; their results are presented below. 

In the case when chemical potential is fixed and placed within the electron subbands 
(as one can see comparing figures ^ja and §3 for /i = —0.4) taking into account of 
fluctuations does not change qualitatively the results obtained previously within the self- 
consistent GRPA [|13| (except for some special cases presented below) . The quantitative 
changes due to fluctuations are important in the region of the critical point leading 
to signiflcant lowering of the critical point temperature (T=0.018 for the self-consistent 
GRPA in flgure ^ and T=0.0145 for the Onyszkiewicz type approximation in flgure ^d). 

Temperature behaviour of the pseudospin mean value and r.m.s. fluctuation 
parameter is presented in flgure ^for the flxed value of the longitudinal fleld h = 0.204. 
At increasing temperature, the pseudospin mean value and r.m.s. parameter jump from 
the branch of the low temperature phase to the one of the high temperature phase in 
the phase transition point T = 0.0109. 

A certain qualitative change takes place (flgure |^) when the chemical potential is 
placed near the center of electron subband (yU = —g/2 = —0.5). Comparing to the self- 
consistent GRPA (flgure |^a) or the DMFA (flgure ^d), a change of slope of the phase 
coexistence line is observed. The vertical line on the — h phase diagram (flgure |^) 
twice crosses the phase coexistence curve and hence there exists the possibility of the 
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Figure 4. Temperature dependences of (a) pseudospin mean value and (b) r.m.s. 
fluctuation parameter {g = 1, <fc=o = 0.2, /i = —0.4, h = 0.204). Thick solid lines 
correspond to thermodynamically stable states, other lines denote metastable and 
unstable ones. 
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Figure 5. Tc — ft. phase diagrams {g — 1, tk=o — 0.2, /i = —0.5). 



two sequential first order phase transitions (reentrance) with change of temperature. 

For the different values of chemical potential, a slope of the phase coexistence curve 
can vary (/i = —0.5 in figure ^ and fi = —0.4 in figure |^). Within the region, where 
change of slope occurs (/x ~ —0.6 ±0.005 and fi ^ —0.425 ±0.005), a possibility of three 
sequential phase transitions of the first order at change of temperature is observed (phase 
diagram in figure ^a). Near the region of the critical point, a qualitative change (change 
of topological type of the phase coexistence curve and appearance of the triple point are 
observed) due to gaussian fluctuations becomes significant when chemical potential has 
a value within the mentioned range (figure Pd). 

The field dependences of pseudospin mean value and r.m.s. fluctuation parameter 
in the temperature vicinity of this region (near the triple point) are presented in figure ^ 
and in figure |^ ((a) - above the triple point, (b) - below the triple point). The phase 
transition points (denoted as a, b and c points in figures) correspond to the crossing 
points of different branches of Q{h) (figure^. 
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Figure 6. Tc — h phase diagrams: (a) self-consistent GRPA, (b) gaussian fluctuation 
method {g = 1, tfc=o — 0.2, fi = —0.6). 



In the figures presented above the case when chemical potential is placed in the 
lower subband is presented. If the chemical potential is placed in the upper subband 
our results transform according to the mentioned above electron-hole symmetry of the 
initial Hamiltonian. 

5. Conclusion 

Self-consistent method, taking into account corrections due to gaussian fluctuations of 
effective field in the self-consistent GRPA scheme, and its simplified version (analogous 



to the approximation proposed by Onyszkiewicz [|T^, ^ for spin models) are apphed 
for investigation of pseudospin-electron model. 

Calculations of the thermodynamic functions have shown that Onyszkiewicz type 
approach (taking into account the class of diagrams restricted as against the usual 
gaussian fluctuation approach) does not change qualitatively any of the results obtained 
within the gaussian fluctuation approximation. A comparison with mean field type 
approximations (self-consistent GRPA, DMFA) shows that taking into account of 
fluctuations is essential in the region of the critical point and can lead to the qualitative 
changes of behavior of thermodynamical functions and the shape of corresponding phase 
diagrams for some special value of chemical potential. The presented results demonstrate 
that the quantitative changes due to fluctuations lower the critical point temperature 
and shift the corresponding value of longitudinal field when chemical potential is 
placed in electron bands. The lowest value of critical temperatures correspond to the 
Onyszkiewicz type approximation. 

Preliminary analysis of temperature behaviour of the pseudospin correlator {S^S^)q 
(p^) (with fixed r.m.s. fluctuation parameter) shows that the high temperature phase 
become unstable with respect to the fluctuations with q when chemical potential 
is placed between the electron subbands. The maximal temperature of instability is 
achieved for q = (tt, vr) and indicates a possibility of phase transition into a modulated 
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Figure 7. Field dependences of pseudospin mean value: (a) above the triple point 
(T = 0.01370), (b) below the triple point (T = 0.01316) {g = 1, ife=o = 0.2, /i = -0.6). 




Figure 8. Field dependences of r.m.s. fluctuation parameter: (a) above the triple 
point (r = 0.01370), (b) below the triple point (T = 0.01316) {g = 1, tk=o = 0.2, 
/i = -0.6). 




Figure 9. Field dependences of grand canonical potential: (a) above the triple point 
(T = 0.01370), (b) below the triple point (T = 0.01316) {g = 1, tk=o = 0.2, /x = -0.6). 
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Figure 10. Temperature dependences of the pseudospin mean value: (a) h = 0.01467, 
(b) h = 0.01495. 



(chess-board) phase. This resuh confirms the previously obtained within the framework 
of the self-consistent GRPA one [|r3|, 0, but taking into account of fluctuations lowers 
value of temperature in which the instability occurs. 

In this paper we investigated the possible phase transitions in PEM within the 
gaussian fluctuation approximation without creation of super structures {q = 0) and all 
the presented phase diagrams concern only this case. Presented in our paper method 
of consideration of gaussian fluctuations of the effective mean field can be used to 
investigation of the influence of fluctuations on thermodynamic properties of modulated 
(chess-board) phase (like it was done in work within the framework of the self- 
consistent GRPA). This issue will be the subject of the further investigation. 
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